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Mapping brain glucose uptake with chemical exchange- sensitive 
spin-lock magnetic resonance imaging 

Tao Jin\ Hunter Mehrens^ Kristy S Hendrich^ and Seong-GI Kinn^'^'^ 

Uptake of administered D-glucose (Glc) or 2-deoxy-D-glucose (2DG) has been indirectly mapped through the chemical exchange 
(CE) between glucose hydroxyl and water protons using CE-dependent saturation transfer (glucoCEST) magnetic resonance 
imaging (MRI). We propose an alternative technique — on-resonance CE-sensitive spin-lock (CESL) MRI — to enhance responses to 
glucose changes. Phantom data and simulations suggest higher sensitivity for this 'glucoCESL' technique (versus glucoCEST) in the 
intermediate CE regime relevant to glucose. Simulations of CESL signals also show insensitivity to 8o-fluctuations. Several findings 
are apparent from in vivo glucoCESL studies of rat brain at 9.4 Tesia with intravenous injections. First, dose-dependent responses 
are nearly linearly for 0.25-, 0.5-, and 1-g/kg Glc administration (obtained with 12-second temporal resolution), with changes 
robustly detected for all doses. Second, responses at a matched dose of 1 g/kg are much larger and persist for a longer duration for 
2DG versus Glc administration, and are minimal for mannitol as an osmolality control. And third, with similar increases in steady- 
state blood glucose levels, glucoCESL responses are ^2.2 times higher for 2DG versus Glc, consistent with their different metabolic 
properties. Overall, we show that glucoCESL MRI could be a highly sensitive and quantifiable tool for glucose transport and 
metabolism studies. 
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INTRODUCTION 

Glucose is the body's fuel; its uptake and metabolism are sensitive 
biomarkers of cellular dysfunction in afflictions such as cancer, 
stroke, Alzheimer's, and psychiatric diseases. Direct detection of 
glucose uptake in vivo has been established in human and animal 
studies by ^^C and nuclear magnetic resonance spectroscopy,^""^ 
but low glucose concentrations severely limit sensitivity. Clinical 
positron emission tomography scans with radioisotope-labeled 
deoxyglucose detect altered metabolism, but radiation exposure 
limits the scan frequency and excludes certain patient groups. In 
addition, this procedure is expensive and has low spatial and 
temporal resolutions. Thus, development of magnetic resonance 
imaging (MRI) techniques sensitive to levels of nonradioactive 
glucose and its analogs is of crucial importance. Throughout this 
paper, we refer to glucose and its analogs simply by the term 
'glucose'. 

Many key biomolecules can be indirectly detected via their 
chemical exchange (CE) with water. One MRI approach uses the CE 
saturation transfer (CEST) technique,^"^ where low-power, long- 
duration (several seconds) radiofrequency (RF) irradiation at the 
resonance frequency of labile protons (e.g., the hydroxyl protons 
in glucose) provides a magnetic label that is transferred to water 
protons by CE. This saturation transfer attenuates the bulk water 
signal, improving the detectability to labile protons. Recently, 
animal CEST MRI studies administering nonradioactive D-glucose 
(Glc) and 2-deoxy-D-glucose (2DG) have shown promising results 
for tumor detection and cell-line differentiation,^ and for 
assessment of glucose uptake and metabolism in normal brain. 



These encouraging studies indicate that 'glucoCEST' has a wide 
range of applications; however, as will be shown below, technical 
challenges include relatively low sensitivity, susceptibility to Bq 
inhomogeneity and fluctuation, low temporal resolution, and 
difficulty in quantification. A main concern is that optimal CE 
contrast is obtained by CEST only when the exchange rate (k) 
between protons in water and those in the hydroxyl groups of 
glucose is much slower than their resonance frequency separation 
(df'^"^ (i.e., k/^«^); however, this condition is not satisfied 
under physiologic conditions where exchange between water and 
glucose hydroxyl protons is rapid, and their resonance frequencies 
are close.^^ Therefore, high glucose doses are required for CEST 
detectability in normal brain. 

An alternative technique to indirectly map the CE of biomole- 
cules with water utilizes on-resonance spin-lock (SL) MRI. This 
CE-sensitive SL (CESL) technique offers promise as a tool to 
enhance sensitivity to faster exchange processes (versus 
CEST).^^'^^ We recently reported that CESL contrast is optimal for 
the intermediate CE regimes relevant to glucose (i.e., k/3 ~ 1), and 
showed high glucose detection sensitivity in phantom studies.^^ In 
addition, feline and human applications of SL functional MRI have 
shown high-temporal resolution, which can improve the 
detectability to glucose changes.^^'^° On-resonance SL signals 
are usually quantified by Tip, the spin-lattice relaxation time in the 
rotating frame. Although Tip has contributions from all CE and 
non-CE relaxation pathways, and thus lacks specificity to any 
particular type of labile protons, recent in vivo studies have shown 
that tissue Tip is insensitive to variations in blood flow, volume. 



^Department of Radiology, University of Pittsburgh, Pittsburgh, Pennsylvania, USA; ^Department of Physics, University of Pittsburgh, Pittsburgh, Pennsylvania, USA; ^Center for 
Neuroscience Imaging Research, Institute for Basic Science, SKKU, Suwon, Korea and "^Department of Biological Sciences, SKKU, Suwon, Korea. Correspondence: Dr T Jin, 
Department of Radiology, University of Pittsburgh, 3025 E Carson Street, Pittsburgh, PA 15203, USA. 
E-mail: taj6@pitt.edu 

This work is supported by NIH grants EB003324, P30-NS076405, and T90-DA022761. 

Received 30 January 2014; revised 8 April 2014; accepted 10 May 2014; published online 28 May 2014 



Mapping glucose uptake in brain 
T Jin et ol 



and intravascular susceptibility;^ ^'^^ this lack of chennical specificity 
(versus CEST) is alleviated in applications with glucose adnnini- 
stration (referred to as 'glucoCESL'), where the the signal response 
is expected to be donninated by elevated concentrations of 
glucose and its nnetabolic products. 

In the current work, we nnethodically exannined the feasibility of 
mapping brain glucose changes with CESL after intravenous 
injection of Glc or 2DG. Assessnnent included 9.4 T (Tesia) 
experinnents and sinnulations incorporating SL frequency variation. 
Phantonn experinnents of bionnolecular solutions characterized 
CESL for hydroxyl groups with different k values, and for varying 
pH, concentration, and water longitudinal (T^) and transverse (T2) 
relaxation tinnes. Connputer sinnulations of CESL versus CEST 
connpared sensitivity features and investigated the dependence 
on Bq shifts. Three separate groups of in vivo rat-brain glucoCESL 
studies with high-tennporal resolution were perfornned with 
intravenous injections. First, the dose dependence and detection 
threshold with our settings was deternnined for Glc. Next, 
glucoCESL responses were connpared for nnatched doses of Glc, 
2DG, and nnannitol, and with nneasurennent of arterial blood Glc 
and 2DG levels. The latter studies served to characterize dynannic 
properties and to assess A/?ip contributions fronn changes in tissue 
water content after intravenous injection of hypertonic solutions; 
nnannitol served as a control, since it induces large elevations in 
plasnna osnnolality, has a long half-life in blood (~60 nninutes in 
rats^^), and does not cross the blood-brain barrier.^^ And lastly, 
glucoCESL responses were obtained for varied SL power levels at 
steady-state blood glucose levels to further exannine the origin 
and quantification of glucoCESL responses. 

MATERIALS AND METHODS 

Theoretical Background 

Although CEST and CESL each incorporate an RF preparation period 
to encode the exchange processes between water and labile protons, 
their labeling approaches are different. In the on-resonance CESL 
technique, water proton spins are excited (i.e., magnetization is flipped 
to the transverse plane) and then SL irradiation is applied at that same 
water resonance frequency. During the SL period, water protons exchange 
with labile protons (e.g., the hydroxyl groups of glucose). Water proton 
magnetization is then returned to the longitudinal plane, where the CE- 
reduced magnetization can be detected. The SL frequency in CESL is 
usually a few hundred Hertz (Hz) to effectively 'lock' the water 
magnetization in the presence of Bq inhomogeneity, and this requires 
higher power (equivalent to an RF magnetic field strength B-i of ~10^T) 
versus CEST studies, where B-i is usually <2iuJ. However, the 
on-resonance SL duration in CESL is usually comparable to water T2 (on 
the order of tens of milliseconds for in vivo studies), and this duration is 
much shorter than required for the off-resonance saturation of CEST, which 
should be comparable to water Ti (on the order of seconds). Thus, the 
specific absorption rate for CESL is often similar to that of CEST, and 
has already been applied to several human studies.^^'^"^ 

Signal intensity in CESL is determined by /?ip (= ^/T-^p), the spin-lattice 
relaxation rate in the rotating frame, and /?Tp is obtained by varying the SL 
duration (TSL), and fitting the data to 5 = So • exp(-/?Tp • TSL), where 5 and 
So are signal intensities with and without SL irradiation, respectively. 
Assuming a two-pool exchange process where the relative population of 

labile to water protons, p, is much smaller than 1, /?ip can be expressed 

3318,25 

f^^p = /?ex + /?2,0, [1] 

where Rex is the exchange-mediated relaxation rate, and /?2,o is the 
transverse relaxation rate of water protons without CE contributions. The 
Rex term fully characterizes the CE process and is proportional to labile 
proton concentration as defined by^^'^^ 

Re.= ,'-'-f , [2] 

where co-] is the SL frequency equal to yB-^/ln (units of Hz), and y is the 
gyromagnetic ratio. Determination of R^p at multiple coi levels yields a 
'dispersion curve', where /?2,o can then be approximated as the /?ip value at 
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high CO-], and CE-related values p, /c (units of per second), and d (units of 
rad/second) can be determined by fitting the dispersion data to Equations 
[1] and [2].^^ Equation [1] quantitatively describes both CESL and CEST data 
for glucose and other biomolecules, as confirmed by our recent phantom 
studies.^^ 

In glucoCESL applications, A/?ip with administration of Glc or 2DG should 
be dominated by ARex, and relate to the change in the ratio of hydroxyl to 
water protons (Ap), as expressed by 

AR^p = A/?ex + A/?2,o ^ ARex = ,2 ^^. o l ' [^l 



Sinnulations 

Characterizations of CESL and CEST were performed by computer 
simulations with Bloch-McConnell equations as follows. 

(1) Sensitivity to different exchange rates was compared for CESL versus 
CEST by simulating the maximum CE contrast achievable by each of these 
approaches for various l</d values, for water /?i = 1 /Ti = 0.5 per second, and 
water /?2,o values of 5 and 20 per second. For comparison with in vivo 
glucoCESL data, we also assumed a chemical shift, d, between water and 
labile protons of 3,770 rad/second, equivalent to +1.5ppm relative to 
water at 9.4 T, and a labile proton concentration of lOmmol/L, with 
water /?i =0.5 per second and water /?2,o = 20 per second. Contrast for 
CESL was calculated as the difference (AS) between Tip-weighted signals 
(5) with CE (i.e., Pt^O) and without CE (i.e., p = 0), normalized by So. 
Contrast for CEST was obtained from the asymmetric magnetization 
transfer (MT) ratio, which is the difference (AS) between signals with 
irradiation frequency offsets, Q, at the reference frequency (Q= -d) and 
at the hydroxyl proton frequency (Q= +d), normalized by So- Maximum 
contrast at each k value was simulated for CESL at optimal coi and TSL 
values, and for CEST at optimal a>i for steady-state conditions (saturation 
duration^ 00). 

(2) The effect of errors on CESL versus CEST measurements owing to 
variations in water resonance frequency was simulated for Bq shifts ranging 
from 0 to 20 Hz, assuming water /?i=0.5 per second, water /?2,o = 20 
per second, and no CE (p = 0). The CEST asymmetric MT ratio errors were 
calculated at Q = 200, 400, and 600 Hz with fii = 1.6/^T and 4-second 
saturation duration (typical for glucoCEST studies^°'"'^^). For comparison, 
the errors in CESL contrast (AS/So) and CESL R^p owing to Bq shifts were 
simulated for co^ =300, 400, and 500 Hz. Since Bq shifts cause CESL signal 
oscillations, the magnitude of signal oscillations was first calculated as 
AS = |STSL(with Bq shift) -STSL(without Bq shift)| for each of 100 equally 
spaced TSL values ranging between 45 and 55 milliseconds (typical in vivo 
tissue 7"ip values at 9.4T^^) and then averaged. 

Chennical Exchange-Sensitive Spin-Lock Magnetic Resonance 
Innaging Experiments — General 

All MRI experiments were performed on a 9.4-T MR system (Agilent 
Technologies, Santa Clara, CA, USA). For phantom experiments, a 3.8-cm ID 
quadrature volume coil (Rapid Biomedical, OH, USA) provided RF 
transmission and detection; on-resonance SL irradiation of TSL duration 
at the desired coi level was sandwiched between a 90° hard excitation 
pulse and a 90° hard flip-back pulse. For the in vivo experiments of 
paradigms 1 and 2 (see below), RF transmission was provided by a volume 
coil (6.4 cm diameter) and received by a surface coil (2.2 cm diameter; coil 
combination from Nova Medical, MA, USA). For the in vivo experiments of 
paradigm 3, RF was transmitted and received by a home-built surface coil 
(1.6cm diameter), since these studies required higher coq values. In all 
in vivo studies, 50 milliseconds of SL irradiation at the desired cdt level 
was sandwiched between an adiabatic half-passage pulse and an 
inverse adiabatic half-passage pulse (each of 1.5 millisecond duration).^^ 
Immediately after SL preparation, MR images were acquired by single-shot 
echo planar imaging (EPI). 

Chemical Exchange-Sensitive Spin-Lock Phantom Magnetic 
Resonance Imaging Studies 

Biochemicals dissolved in phosphate-buffered saline were placed in 
individual vials, and data from multivial phantoms were acquired at 37°C 
to characterize various CESL properties, (i) Effects of hydroxyl groups from 
different molecules were studied with 30 mmol/L solutions of Glc, 2DG, 
glycogen, and myo-inositol at pH = 7.0. (ii) Dependence on pH was 
determined with 20 mmol/L Glc at pH = 6.8, 7.0, 7.2, and 7.4. (iii) Water 
and T2 effects on quantification were examined with 5 and 20 mmol/L 



© 2014 ISCBFM 



Journal of Cerebral Blood Flow & Metabolism (2014), 1402-1410 



Mapping glucose uptake in brain 
T Jin et al 



1404 

solutions of Glc at pH = 7.0, both with and without addition of 0.1 mmol/L 
MnCb. (iv) Linear dependence of /?ip on glucose concentration was 
validated with 5, 10, 15, 20, 30, and 40 mmol/L solutions of Glc and 2DG at 
pH = 7.0, with addition of 0.1 5 mmol/L N\nC\2 to reduce T2 to values closer 
to tissue, and also to speed acquisition. For all phantoms, dispersion data 
were acquired with spin-echo EPI for 11 co-i values between 125 and 
4,000 Hz, with 64x64 matrix size, 4x4cm^ field of view, 5mm slice 
thickness, echo time = 25 milliseconds, and repetition time = 15 seconds. 



Chemical Exchange-Sensitive Spin-Lock In Vivo Magnetic 
Resonance Imaging Studies 

Animal preparation. With approval by the Institutional Animal Care and 
Use Committee at the University of Pittsburgh, 32 male Sprague-Dawley 
rats were studied (26 for MRI experiments, and 6 for bench-top 
measurements of blood glucose). Animals were anesthetized with 
isoflurane (5% for induction and 2% for surgery) in a mixture of O2 and 
air, with the O2 concentration maintained at 30%. The right femoral vein 
was catheterized to deliver maintenance fluid and Glc (or 2DG or 
mannitol). The right femoral artery was catheterized to ensure arterial 
blood pressure was maintained within a normal physiologic range. After 
surgery, the isoflurane level was reduced to ~ 1 .5%, and maintained at this 
level. End-tidal CO2 levels were monitored and kept within 3.0 to 4.0%. 
Rectal temperature was maintained at 37.2±0.5°C using a feedback- 
controlled heating pad. 

Animal magnetic resonance imaging experiments. GlucoCESL rat-brain 
studies were performed for three separate paradigms with intravenous 
injection of Glc, 2DG, or mannitol. General EPI parameters were matrix 
size = 64x64, field of view = 3.2x3.2 or 2.56 x 2.56 cm^ 2 mm slice 
thickness, repetition time = 3 seconds, echo time = 30 milliseconds for 
single- or double-spin-echo EPI, and echo time = 12.5 milliseconds for 
gradient-echo EPI. Determinations of /?Tp time courses in each run were 
made from interleaved acquisitions of images with TSL = 0 and with 
TSL = 50 milliseconds at the appropriate coi values; when TSL was 0, delays 
after data acquisition were increased by 50 milliseconds to maintain 
constant repetition time. Solutions of Glc, 2DG, and mannitol were 
prepared at 20% weight by volume, and each dose was injected over a 
1- to 5-minute period, except for the additional infusion of Glc in paradigm 3. 
Specific details are as follows: 

Paradigm 1: The detection threshold of glucoCEST MRI with our settings 
was determined by administering consecutive doses of 0.25, 0.5, and 1 .Og/kg 
body weight Glc (n = 5 rats) at 40, 100, and 180 minutes, respectively, after 
initiation of scanning. The 60- and 80-minute intervals between doses 
allowed blood glucose to recover to baseline levels (as shown below in 
Figure 5D for 1 g/kg injection). Maps of R^p (and A/?ip) were generated 
from spin-echo EPI images without SL preparation (i.e., TSL = 0, single 
acquisition) and SL-weighted images with cl>i=500Hz (TSL = 50 milli- 
seconds, 3 averages), yielding a temporal resolution of 12 seconds. 

Paradigm 2: The dynamic properties of glucoCESL MRI were assessed for 
Glc, 2DG, and mannitol administration (n = 4 rats each) for matched doses 
of 1 g/kg body weight, both to compare sensitivity, and to characterize the 
source of AR-ip changes. Data were acquired for 100 minutes after Glc or 
2DG administration; responses after mannitol injection were small, so 
acquisition in those studies was terminated at 60 minutes. Maps of /?ip 
(and A/?ip) were generated as in paradigm 1. In separate bench-top 
studies, blood glucose levels were sampled every 10 minutes with 
administration of Glc (n = 3 rats) or 2DG (n = 3 rats), and time courses were 
compared with those of AR-ip. 

Paradigm 3: The origin of glucoCESL contrast was further assessed by 
acquiring time-dependent data for A/?ip at varied co-i values in Glc and in 
2DG studies with data acquired at steady state. If the origin of signal 
change is dominated by CE effects, then A/?ip should decrease with co-i, as 
expected in Equation [3]. A second goal was to compare steady-state 
responses to 2DG and Glc administration, where the correlations between 
blood and brain glucose levels for the latter have been established.^ As 
blood glucose levels decrease faster after Glc versus 2DG injection, steady- 
state levels were attained in Glc studies (n = 4 rats), by injecting 0.3 g/kg 
body weight Glc over a 1 -minute period, followed by constant infusion of 
2.0 g/kg per hour Glc for 1 hour; four double-spin-echo EPI Tip-weighted 
images were acquired in an interleaved manner, one without and three 
with SL weighting (coi=500, 1,000, and 2,000 Hz). Blood glucose levels 
were sampled shortly before and ~60 minutes after onset of the first Glc 
injection. In 2DG studies (n = 5 rats), a single injection of 1 g/kg 2DG was 
administered; three Tip-weighted images were acquired in an interleaved 



manner either by double-spin-echo or gradient-echo EPI, one without SL 
preparation and two with SL weighting of coq =500 and 2,000 Hz (results 
from both EPI techniques were similar, so they were combined). To ensure 
that coq was similar in all studies and across the region of interest (ROI), a 
midcortical ROI in the dorsal area close to the surface coil was first selected 
from anatomic images (see orange pixels of Figure 5C inset image for an 
example), and was adjusted so that its spatial average over the ROI 
reached the target value.^^ 

Data Analysis 

Data were analyzed with in-house Matlab programs and STIMULATE 
software.^^ For phantoms studies, quantitative analyses were performed 
for large circular ROIs defined within each sample vial. Values of d, k, and 
/?2,o for each solution were first obtained by fitting the /?Tp dispersion data 
to Equations [1] and [2]. Then /?ex values were determined at 0*1 =500 Hz 
and at coi = 2,000 Hz, and the ratio between them was calculated. For 
fitting to Equation [2], the relative population (p) is the number of labile 
protons per biomolecule (N) times its concentration divided by IIOmol/L 
(proton concentration in the water pool), where N is 5 for Glc, 4 for 2DG, 3 
for glycogen, and 6 for myo-inositol. Although different hydroxyl protons 
within each molecule have their own chemical shifts and exchange 
rates^'^^'^^ (e.g., Glc hydroxyl protons shows distinct resonances at ~1.2, 
2.2, and 2.8 ppm relative to water^), it is the average 3, k, and /?ex values 
from all hydroxyl resonances that were determined here. 

For in vivo experiments where images with varied co-i and TSL values 
were acquired in an interleaved manner, images with same coq and 
TSL values were extracted from each run. Time series of R^p images 
were calculated by mono-exponential fitting of signal intensities from 
TSL = 0 milliseconds and TSL = 50 milliseconds, as R^p = \n (Stsl = o/ 
Stsl=5o miiiiseconds)/50 milNseconds. Statistical maps corresponding to /?ip 
changes resulting from Glc or 2DG injection were calculated by performing 
Student's f-tests on a pixel-by-pixel basis, where a threshold of P<0.05 
(uncorrected for multiple comparisons) was chosen and a minimum cluster 
size of 3 pixels was applied. Baseline data were averaged from the period 
spanning 10 to 0 minutes before start of injection. Data for glucose-induced 
changes were from periods spanning 10 to 25 minutes post injection for 
paradigms 1 and 2, and 50 to 70 minutes post injection for paradigm 3. 
Quantitative analyses were performed for ROIs encompassing all brain pixels 
(paradigm 1) and from small midcortical ROIs (paradigms 2 and 3), as 
determined from anatomic images. Time-dependent series of A/?ip maps are 
presented as averages over 10-minute intervals before subtraction from the 
baseline average, while A/?ip time courses from within ROIs are presented as 
averages over 1 -minute intervals and are shown as mean±s.e.m. 

RESULTS 

Chemical Exchange-Sensitive Spin-Lock Magnetic Resonance 
Imaging Studies in Phantoms 

Figure 1 shows dispersion data characterizing CESL properties in 
phantoms. In Figure 1A, 30 mmol/L solutions of Glc, 2DG, myo- 
inositol, and glycogen all show high /?ip values when coi ^500 Hz, 
but this CE effect is attenuated as coi increases, and /?ip is much 
slower when coi ^ 2,000 Hz; this coi dependence on CE is an 
important property of CESL. The arrows in the dispersion curve 
approximately indicate the half-width at half-maximum value of 
coi, where the condition co? = /c^ + ^^ is satisfied. Assuming 
similar S, the higher half-width at half-maximum values for Glc and 
2DG suggest that k is faster than for glycogen and myo-inositol. 
Figure IB shows /?ip dispersion data for 20 mmol/L Glc solutions 
spanning the range of physiologic pH values (6.8 to 7.4); the half- 
width at half-maximum values of coi increase with pH. Character- 
istics of Figure IB dispersion curves can be appreciated from 
Equation [2], given that 3 is constant for all four solutions 
(pH-independent), whereas k increases with pH. Namely, at very 
small coi values, the largest /?ip is seen for the pH 7.0 solution, 
meaning that 3 is closest to the k value of that solution (i.e., in the 
intermediate CE regime), whereas at very large coi, R^p increases 
with /c.^^'^^ Note that at co^ -500 Hz (vertical dashed line), /?ip, and 
therefore Rex, are relatively insensitive to this physiologic pH 
range, with a variation of <8%. However, when coi ^2,000 Hz, /?ip 
and /?ex are much more sensitive to pH. From Equation [2], the 
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Figure 1. Chennical exchange-sensitive spin-lock nnagnetic resonance innaging data at 9.4 T (Tesia) fronn biochennical solutions prepared in 
phosphate-buffered saline, and nneasured at 37°C. (A-C) Dispersion curves (rotating-franne spin-lattice relaxation rates, R^p, versus spin-lock 
frequencies, coi) are displayed for solutions with (A) D-glucose (Glc), 2-deoxy-D-glucose (2DG), glycogen (Gly), and nnyo-inositol (Ins), all at 
30nnnnol/L concentration at pH = 7.0; (B) 20nnnnol/L Glc at four pH values; and (C) 5 and 20nnnnol/L Glc at pH = 7.0 without and with MnC\2. 
Color-nnatched arrows in panel A indicate the approxinnate half-width at half-nnaxinnunn (HWHM) values of coi for each solution, suggesting 
that the chennical exchange rates (/c) between protons in water and those in the hydroxyl groups are faster for Glc and 2DG versus Gly and Ins 
(see text for rationale). The dispersion curve in panel B is steepest for sannples at pH = 6.8 and 7.0, with HWHM values for coi increasing with 
pH value; the /?ip values are sinnilar for all pH values at coi =500 Hz (vertical dashed line). The nearly identical /?ip span between arrows at 
coi =500 Hz in panel C indicates that Glc concentration-dependent R^p differences are independent of and T2 relaxation. (D) As expected 
fronn Equation [2], /?ip is shown to be linearly proportional to glucose concentration in Glc and 2DG solutions (at pH = 7.0 with 0.15nnnnol/L 
l\/lnCl2; acquired with coi = 500 Hz); the slope of the linear fit is 0.066 per second per nnnnol/L for Glc and 0.050 per second per nnnnol/L for 2DG. 



Table 1. Values^ derived by fitting phantom R^p 


dispersion data to Equations [1] and [2] 
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groups; S, resonance frequency separation between protons in water and those in the hydroxyl groups; Rex, 


exchange-mediated relaxation term, ^Results are 
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ratio between /?ex acquired at 500 Hz versus 2,000 Hz should 
therefore decrease nnonotonically with increasing k. Information on 
k (and pH) changes can thus be inferred fronn this ratio of /?ex values. 
Table 1 shows results of fitting /?ip dispersion data to yield values for 
k, 3, and the ratio of /?ex values at coi = 500 versus 2,000 Hz. 

Dispersion data were also acquired for Glc concentrations of 5 
and 20nnnnol/L, with the water relaxation rates altered by IVlnCl2. 
Figure 1C shows that glucose concentration-dependent differ- 
ences in /?ip are independent of the l\/lnCl2-induced differences in 
water Ti and Tj, as seen by the sinnilar spans indicated by the blue 
arrows. Therefore, AR^p fronn Equation [3] can be a quantitative 
index to determine changes in glucose concentration. The linear 
dependence of /?ip on glucose concentration was further 
confirmed by studies of Glc and 2DG solutions with pH 7.0 at a 
typical co-i value of 500 Hz (Figure 1 D). A linear fit to concentration 
yields A/?ip (and A/?ex) of ~ 0.066 per second per mmol/L Glc 
(r = 0.998) or 0.050 per second per mmol/L 2DG (r = 0.999). 



Simulations of Chemical Exchange Effects in Chemical Exchange 
Saturation Transfer and Chemical Exchange-Sensitive Spin-Lock 
with Bloch-McConnell Equations 

Values for maximum contrast with CEST and CESL were calculated 
for different k/S values. Contrast with CEST is optimal for slow CE 
regimes {k/S«^), but very low for fast CE regimes {k/3>^). 
However, CESL contrast is optimal for intermediate CE regimes 



{k/S ~ 1 ), and is higher than CEST when k/3>^. Contrast with both 
CESL and CEST decreases as /?2,o increases. For a more specific 
example relevant to glucose in vivo, the /?2,o = 20 per second data 
from the dotted region of Figure 2A were replotted as a percen- 
tage change in Figure 2B, where both Glc and 2DG are repre- 
sented at 9.4 T by the blue vertical line {k/S = 1.29, from Table 1), 
and at 3 T by the green vertical line {k/S = 4.05). These results 
suggest that CESL MR! is a good choice for glucose mapping at 
9.4 T, and an even better choice at lower fields where k/S is larger. 

Computed errors owing to water resonance frequency varia- 
tions in CEST and CESL experiments appear in Figure 3. Simula- 
tions for typical in vivo CEST saturation (i.e., 81 = 1.6 fiJ applied for 
4 seconds) show that the asymmetric MT ratio signal at an 
irradiation offset (Q) of 400 Hz has an error of 0.5% for the 2-Hz Bq 
shift often encountered in brain studies,^^'^^ 1.2% for a 5-Hz shift, 
and 5% for a 20-Hz shift (Figure 3A). Simulations for typical in vivo 
CESL studies (i.e., TSL ~50 milliseconds) show that errors from Bq 
shifts are greatly reduced (versus CEST) for coi ^300 Hz, with negli- 
gible effects on either contrast (A5/5o; Figure 3B) or on calculated 
/?ip values (Figure 3C), even for Bq shifts as large as 20 Hz. 

Dose-Dependence and Detection-Threshold Evaluation of 
glucoCESL with Glc in Rat Brain 

In vivo glucoCESL data from consecutive Glc injections adminis- 
tered at three different doses appear in Figure 4 (paradigm 1). 
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Figure 2. Sinnulations by Bloch-McConnell equations showing 
enhanced sensitivity to glucose detection for chennical exchange- 
sensitive spin lock (CESL) versus chennical exchange saturation 
transfer (CEST). Values shown are for the nnaxinnunn contrast 
achievable. (A) Nornnalized contrast (AS/SoAp • S)) with CEST highest 
in the slow-exchange reginnes (i.e., when /c/(resonance frequency 
separation, 3) where it outperfornns CESL. Nornnalized contrast 
with CESL highest for internnediate-exchange reginnes {k/S ~ 1), and 
CESL outperfornns CEST in the fast-exchange reginnes {k/3>^). 
Contrast in both CEST and CESL depends on /?2,o (the transverse 
relaxation rate in the absence of exchange effects), but scales 
linearly with 3 and with relative population (i.e., glucose hydroxyl 
protons to water protons, p) when p^3«^, so contrast is 
nornnalized by {p-3). (B) Sinnulations are plotted over a snnaller 
range of k/3 values, for reasonable in vivo values of /?2,o = 20 
per second, 1 0 nnnnol/L of labile protons (P = 0.00091 ), and 3 = 3,770 
rad/second; vertical lines indicating the k/3 values appropriate for 
both D-glucose and 2-deoxy-D-glucose at 9.4 T (Tesia; blue line) and 
at 3T (green line) show that CESL outperfornns CEST at both field 
strengths; however, the sensitivity enhancennent is substantially 
higher at ST. 



Figure 4A shows statistical maps calculated from the A/?ip time 
series for each Glc dose in two separate animals. At the lowest 
dose of 0.25 g/kg Glc, increases in /?ip already appear above the 
statistical threshold for many pixels, but with injection of 0.5 g/kg 
Glc, a widespread increase in /?ip is detected, and at 1 g/kg Glc, the 
increase in /?ip is higher yet, with pixels above the statistical 
threshold spanning most of the brain. Figure 4B shows the 
average of the brain /?ip time courses from all five animals, where 
a dose-dependent /?ip increase is robustly observed, even at 
0.25 g/kg Glc. For individual animals, the brain A/?ip value at the 
time of its highest response increases nearly linearly with Glc dose 
(Figure 4C). 

Dynamic glucoCESL Response Comparisons for Glc, 2DG, and 
Mannitol in Rat Brain 

Time-dependent glucoCESL (Figure 5) shows in vivo postinjection 
responses that are much larger and persist for a longer duration 
for 1 g/kg 2DG versus 1 g/kg Glc (paradigm 2) in serial A/?ip maps 
typifying 0 to 90 minutes post injection (Figures 5A and 5B, 
respectively). These dynamic properties are further shown in A/?ip 
time course averages from all the four animals (Figure 5C). Peak 
A/?ip values were reached in ~20 minutes for both 2DG and Glc, 
but with ~2.5 times higher magnitude for 2DG, which is in good 
agreement with a glucoCEST report.^ ^ Recovery of A/?ip to 
preinjection baseline levels occurred around 60 minutes post-GIc 
injection, but recovery for 2DG still had not occurred at the end of 
the 100-minute window. The higher-intensity and longer-duration 
A/?ip responses to 2DG (versus Glc) can partially be explained by 
the difference in time courses of blood glucose levels after single 
1-g/kg dose injections (Figure 5D). Blood glucose levels peaked 
immediately after injection, with magnitudes similar for 2DG and 
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Figure 3. Simulations by Bloch-McConnell equations showing that small magnetic field {Bq) shifts give large errors in chemical exchange 
saturation transfer (CEST), and very small errors in chemical exchange-sensitive spin lock (CESL). Large Bo-dependent errors appear for (A) 
CEST asymmetric magnetization transfer ratio (MTRasym) in simulations for three radiofrequency offset (Q) values, while negligible 
8o-dependent errors are present for (B) CESL contrast (AS/Sq), and (C) CESL /?ip values in simulations at three coi levels. Note the vastly different 
vertical axes scales in panel A versus panel B. 




Ti me { min u ies ) Glucose dose (g/kg) 

Figure 4. Rat-brain glucoCESL studies at 9.4T showing near-linear contrast for intravenously administered D-glucose (Glc) doses of 0.25, 0.5, 
and 1.0 g/kg, and robust detection for doses ^0.25 g/kg (/V? vivo paradigm 1). (A) The f-maps for each dose for two of the animals show 
highest f-values in the cortex where sensitivity is higher with our surface coil reception. Color scale: f-value. (B) Average of time courses for all 
animals {n = 5, mean±s.e.m.) clearly shows the increase in brain A/?ip with Glc dose. Arrows indicate time of injection. (C) The nearly linear 
dependence of peak brain A/?ip on Glc dose appears for each individual animal. CESL, chemical exchange-sensitive spin lock. 
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Figure 5. Tennporal dynannic properties in rat brain showing larger and longer-duration /?ip responses at 9.4T after the intravenous injection 
of 2-deoxy-D-glucose, 2DG (1 g/kg) versus D-glucose, Glc (1 g/kg) — with only nninor contributions fronn osnnolality changes (nnannitol, 
1 g/kg) — and nnore persistent elevation of blood glucose levels for 2DG versus Glc {in vivo paradignn 2). Representative tinne-resolved nnaps 
show A/?ip before and 0 to 90 nninutes after the injection of Glc (A) versus 2DG (B), where tinnes are 0 to 10 and 80 to 90 nninutes for the first 
and final postinjection nnaps, respectively; note the difference in A/?ip gray scale ranges under each series. High-tennporal resolution A/?ip tinne 
courses (C) of Glc, 2DG, and nnannitol {n = 4 each, nnean ± s.e.nn.) are shown for nnidcortical regions as typified by orange pixels in the inset 
innage; tinne courses fronn all brain pixels are qualitatively sinnilar to those of the nnidcortical regions (not shown). The hypertonic nnannitol 
injection serves as an osnnolality control to investigate contributions to A/?ip owing to any changes in tissue water content. (D) Tinne courses of 
blood glucose changes in bench-top studies with injection of 1-g/kg Glc or 1-g/kg 2DG {n = 3 each, nnean ± s.e.nn.) show dynannic 
characteristics that differ fronn A/?ip (connpare to panel C) at initial tinne points (<20 nninutes), but are sinnilar at later tinne points; gray bar 
indicates injection tinne. CESL, chennical exchange-sensitive spin lock. 



Glc. However, blood glucose levels with Glc decreased rapidly, 
returning to baseline at ~60 nninutes post injection, whereas 
levels with 2DG decreased nnuch slower; after a postinjection 
delay of 10 to 40 nninutes, the ratio of blood glucose concentra- 
tion decrease for Glc versus 2DG is qualitatively similar to that for 
A/?ip (Figure 5C versus 5D). The osnnolality change after 1-g/kg 
mannitol injection in the four animals has a smaller effect (versus 
Glc) on the A/?ip time course (Figure 5C). 



Dependence of In Vivo glucoCESL Contrast on Spin-Lock 
Frequency and Blood Glucose Levels 

Further assessment of contributions to in vivo A/?ip was achieved 
at a steady state (i.e., administration tailored to attain similar 
blood glucose concentrations, Glc versus 2DG) with determination 
of SL-frequency dependence (paradigm 3). Time courses in 
Figure 6 for A/?ip at different coi values represent averages 
from all animals administered Glc {n = 4) or 2DG {n = 5). Arterial 
blood glucose levels sampled at 60 minutes post injection confirm 
that increases for Glc versus 2DG studies were similar: 
160±25mg/dL for Glc during MR! experiments, versus 
144±31 mg/dL for 2DG from separate bench-top measurements 
(see Figure 5D). GlucoCESL data show that A/?ip is much smaller 
for co-i =2,000 versus 500 Hz, both for Glc (Figure 6A) and 2DG 
(Figure 6B), which is expected for responses dominated by CE 
effects, as reflected in Equation [3] and shown in Figure 1. The 
averaged ratio of A/?i p at coi = 500 Hz to A/?i p at coi = 2,000 Hz is 
1.75 for Glc and 2.1 for 2DG. The average steady-state A/?ip value 
at coi=500Hz between 50 and 70 minutes post injection is 



0.30 ± 0.03 per second for Glc, and 0.65 ± 0.03 per second for 2DG, 
respectively. 

DISCUSSION 

Enhancement in Glucose Detection Sensitivity for glucoCESL 
versus glucoCEST 

Previous glucoCEST rat-brain studies at 9.4 T found no detectable 
responses for a dose of 0.5 g/kg Glc, and even at 1 g/kg Glc the 
responses were weak and had diminished after 10 minutes, 
despite the fact that blood glucose levels remained hyperglyce- 
mic.^^ Reliable detection of brain glucose with glucoCEST requires 
that blood glucose levels should be raised above baseline by more 
than 300 mg/dL (- 17 mmol/L).^^ With the near-linear relationship 
between brain and blood glucose levels established by nuclear 
magnetic resonance spectroscopy,^'^ these studies suggested a 
brain glucose elevation of 6 to 7 mmol/L as a detection threshold 
for glucoCEST at 9.4 T. However, our in vivo glucoCESL rat-brain 
studies at 9.4 T show robust detection at doses as low as 0.25 g/kg 
Glc, meaning A/?ip ~0.1 per second is the detection threshold for 
our conditions, which corresponds to ~ 1.5 mmol/L increase in 
brain glucose concentration. We reported in a prior fMRI study 
that with sufficient signal averaging to enhance the sensitivity, 
neural stimulation induced a small CESL A/?ip response of ~0.05 
per second in the cat visual cortex, which may be attributed to a 
submillimolar decrease in glucose.^^ 

This enhancement in glucose detection sensitivity for glucoCESL 
versus glucoCEST may be explained by three major factors — differing 
sensitivity to exchange rates, differing temporal resolution, and 
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Figure 6. Rat-brain glucoCESL contrast at 9.4 T confirmed to be mainly from glucose chemical exchange {in vivo paradigm 3). Blood glucose 
concentrations were sustained at a steady state in (A) by injecting 0.3g/kg D-glucose (Glc) over a 1 -minute duration (blue arrow), followed by 
a constant infusion of Glc at 2g/kg per hour for 1 hour (gray bar). For the 2-deoxy-D-glucose (2DG) studies of (B), a single injection of 1 g/kg 
was given at time = 0. Arterial blood glucose levels at 60 minutes post injection increased to 160±25mg/dL for Glc and 144 ±31 mg/dL for 
2DG. Time courses of A/?i p dependence on coi with intravenous injections of Glc (A, n = 4, mean ± s.e.m.) or 2DG (B, n = 5, mean ± s.e.m.) show 
values within midcortical regions, as typified in the inset image of Figure 5C. For clarity, A/?ip data with coi = 1,000 Hz is not shown, since it falls 
between the data for coi =500 and 2,000 Hz. The A/?ip values are much smaller for coi =2,000 versus 500 Hz, both for Glc and 2DG, which is 
expected for responses mainly due to CE effects. CESL, chemical exchange-sensitive spin lock. 



different effects from Bq shifts. First, k values for Glc were deter- 
mined to be 4,280 to 8,140 per second for a physiologic range of 
pH values (6.8 to 7.4), and this falls between the intermediate and 
fast exchange regimes {k/S>^, Table 1), where CESL contrast is 
expected to exceed that of CEST (see Figure 2B). However, this 
sensitivity advantage of glucoCESL will be reduced if k is much 
slower, such as the case for the lower pH of ischemic tissues. 
Second, temporal resolution in our glucoCESL sensitivity tests was 
12 seconds, while typically for glucoCEST it is ~10 minutes to 
accommodate the long saturation pulse and the need to acquire 
images at multiple RF offsets.^ This much-higher temporal 
resolution of glucoCESL provides a higher statistical power for 
glucose detection as compared with glucoCEST. Third, the effect 
of any Bq shift is much smaller for glucoCESL versus glucoCEST. 
Because the saturation frequency for CEST is quite close to water, 
Bq shift of just a few hertz can cause significant signal errors, 
making an effective pixel-wise Bq correction scheme with 
acquisition of redundant data essential. Unfortunately, 
accurate pixel-wise Bq corrections could be very difficult in the 
presence of the large physiologic or instrumental fluctuations that 
may occur during long CEST acquisitions.^^ 

Sensitivity enhancement of glucoCESL versus glucoCEST is 
expected to be even higher at lower magnetic fields for two 
reasons, (i) Since S is proportional to the magnetic field, the 
increase in k/S at lower clinical fields (such as 3T) as compared 
with ultra-high fields (^9.4T) will consequently make the relative 
sensitivity gain of CESL versus CEST even more prominent 
(Figure 2B). (ii) Although, CESL only have minimal errors origi- 
nating from Bq shifts, these errors in CEST may intensify at lower 
fields where the RF saturation frequency offset is proportionately 
closer to water (see Figure 3A, Q = 200 versus 600 Hz). 

Origin of glucoCESL Signals 

A glucose increase in three separate pools can contribute to 
glucoCESL signals: intravascular blood and extracellular and 
intracellular spaces in tissue. The intravascular blood glucose 
contribution to A/?ip is minimal at 9.4T because of (1) a small 
arterial vascular volume fraction (^1%) and (2) suppression of the 
venous component in EPI acquisitions at relatively long echo 
times, due to short T2 and T2* values.^^ In addition, insight about 
the contribution from intravascular glucose can be obtained by 



comparing the temporal dynamics of blood glucose levels and the 
high-temporal resolution time course of A/?ip. Blood glucose 
peaks immediately after injection, whereas A/?ip increases 
gradually and peaks ~15 to 20 minutes post injection. This 
temporal mismatch indicates that A/?ip should have minimal 
contributions from intravascular glucose. 

A decrease of extravascular water content induced by the 
elevation of plasma osmolality can contribute to a small A/?ip 
term, as evidenced after mannitol injection (Figure 5C). However, 
this contribution is exaggerated with mannitol, since it remains in 
the vasculature with a half-life of ~60 minutes.^^ The osmolality- 
induced A/?ip for glucoCESL is dependent on the difference 
between blood and tissue glucose concentrations. Tissue uptake 
of Glc and 2DG makes their osmolality contributions to A/?ip even 
smaller, since there will be rapid equilibration of blood and tissue 
glucose concentrations. In addition, osmolality effects should be 
further reduced by rapid metabolism of Glc to glucose-6- 
phosphate, and metabolism of 2DG to 2DG-6-phosphate 
(2DG6P) with the latter accumulating in tissue. 

As glucoCESL vascular contributions are minimal, the change in 
/?ip mostly arises from the extracellular and intracellular spaces in 
tissue, where relative contributions depend on glucose concen- 
tration, pool size, k, and coi. Intracellular pH is usually 7.0 to 7.1,^'^^ 
while extracellular pH is ~7.2 to 7.3.^^'^"^ Our phantom data 
suggest that at coi = 500 Hz, glucoCESL sensitivity is similar for this 
range of pH values (Figure IB), and therefore the A/?ip change 
would have similar weighting from both extracellular and 
intracellular spaces and reflect the total glucose change. The ratio 
of in vivo A/?ip at 500 to 2,000 Hz is -40% smaller than for the 
aqueous phantoms at pH = 7.0, suggesting that overall in vivo k 
values are higher (according to Equation [3] and Table 1). It is 
possible that in vivo k values are higher than our phantom at 
matched pH because of a different concentration of exchange 
catalysts. Further studies will be necessary to quantify in vivo k 
values and to determine whether A/?2,o is negligible in Glc and 
2DG results, possibly by fitting in vivo dispersion data acquired 
with more co-i levels to Equations [1] and [2]. 

Quantification of Glucose Concentration Changes 

One critical question is whether a change in Glc concentration can 

be determined with CE-MRI. In both CEST and CESL, the CE 
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sensitivity is coupled with non-CE relaxation nnechanisnns includ- 
ing Ti and T2 effects, and MT fronn sennisolid nnacronnolecules, as 
seen in Figures 2A for varied Tj. Consequently, different and MT 
effects will alter CE signal intensities, even when nnetabolite 
concentrations are the sanne.^^'^^ Therefore, in vivo glucoCEST 
signals cannot be directly connpared with those nneasured in 
aqueous phantonns to deternnine the glucose concentration. 
Unlike CEST and CESL signal intensities, our phantonn results 
indicate that R^p changes in CESL are independent of and T2, 
and provide a quantitative index of glucose concentration 
changes (Figure 1C). In addition, the non-CE MT effect is also 
incorporated into the /?2,o ternn in Equation [1] with nnininnal 
contribution to A/?ip in Equation [3]. Assunning that the /c value is 
the sanne in vivo as for our aqueous phantonns (Figure 1C and 
Table 1), the brain glucose change can be estinnated. In the steady 
state, the A/?ip value of 0.30 per second averaged between 50- 
and 70-nninutes post-GIc injection at coi =500 Hz in paradignn 3 
(Figure 6A) yields an estinnate of 4.5 mnnol/L increase in brain 
glucose concentration, or 3.7/innoL per grann tissue, assunning 
0.83 g water per grann tissue. This is in good agreennent with the 
conversion of our 160±25nng/dL (i.e., 8.9 ± 1.5 nnnnol/L) blood 
glucose levels (obtained at 60 nninutes post injection) to brain 
glucose levels of 4±1 /innoL per grann tissue, using the blood- 
brain glucose relationship.^ 

Although, blood glucose levels were sinnilarly elevated in Glc 
and 2DG studies (160 versus 144nng/dL) when A/?ip reached 
steady state (60-nninute tinne point of Figure 6), the increase in 
brain 2DG concentration is estinnated to be ^13 nnnnol/L fronn 
A/?ip of 0.65 per second, nnuch higher than for Glc. It should be 
noted that glucose-adnninistered CE-MRI techniques detect an 
overall increase in the nunnber of hydroxyl groups, which 
predonninantly arise fronn concentration changes of glucose (Glc 
or 2DG) and its nnetabolic products containing hydroxyl groups, 
such as glucose-6-phosphate and 2DG6P. Given a nnininnal 
contribution fronn intravascular blood, the higher steady-state 
A/?ip values in 1-g/kg 2DG versus 1-g/kg Glc studies (Figure 5) are 
prinnarily because of a larger extravascular (i.e, neural tissue) 
concentration of 2DG and its nnetabolic products (versus Glc). It is 
known that unlike glucose-6-phosphate, 2DG6P accunnulates in 
cells.^^'^^ It was reported that with a 0.5-g/kg 2DG injection, the 
brain concentration of 2DG6P was approxinnately the sanne as that 
of 2DG at 20 minutes post injection and was approxinnately four 
tinnes that of 2DG at 100 nninutes post injection.^^ Thus, with 
sinnilar increases of blood Glc and 2DG concentrations, the nnuch 
higher A/?ip of 2DG nnay be nnainly because of the increased 
concentration of intracellular 2DG6P. This observation and 
interpretation are consistent with the recent glucoCEST results,^ ^ 
where the signal fronn 2DG is nnuch higher than Glc, even when 
the steady-state blood glucose concentration with Glc infusion is 
higher than that of 1-g/kg 2DG injection. 



CONCLUSION 

Recent glucoCEST MRI studies show potential for widespread 
applications in diseases with altered glucose transport and nneta- 
bolisnn, but sensitivity is linnited, especially at clinical nnagnetic 
fields, and quantification is difficult. Our alternative glucoCESL 
approach substantially increases sensitivity and provides an index 
for potential quantitative nnapping of glucose concentration 
changes that can be applied to dynannic studies with subnninute 
tennporal resolution. These advantages nnake glucoCESL MRI a 
viable tool for studies of glucose transport and potentially nneta- 
bolisnn in clinical settings and in functional studies. 
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